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NUMERICAL  EXPERIMENTS  ON  THE 
STRATOSPHERIC-TROPOSPHERIC  DYNAMICAL  INTERACTION 

I.  INTRODUCTION 

Most  tropospheric  weather  prediction  models  impose  an  upper  boun¬ 
dary  condition  of  w  »  3  0  at  some  finite  height  in  the  lower  strato¬ 

sphere.  Lindzen  et  al.  (1968)  pointed  out  that  this  nrigid  lidM  boun¬ 
dary  condition  produces  spurious  reflections  and  introduces  oscillations 
and  standing  waves  whose  amplitudes  can  mask  the  true  solutions.  In  a 
more  recent  study  Kirkwood  and  Derome  (1977)  investigated  the  effect 
of  this  boundary  condition  on  forced  stationary  planetary  waves.  They 
used  realistic  mean  zonal  winds  in  their  quasi-geostrophic  study  rather 
than  a  quiescent,  isothermal  atmosphere  as  Lindzen  et  al.  (1968)  adop¬ 
ted.  The  reflection  and  transmission  of  these  ultralong  waves  are 
strongly  dependent  on  the  structure  of  the  polar  night  jet  centered 
near  the  stratopause  (Charney  and  Drazin,  1961;  Schoeberl  and  Geller, 
1977).  Kirkwood  and  Derome  found  that  the  stationary  planetary  wave 
structure  could  only  be  accurately  modeled  when  the  vertical  resolution 
in  the  stratosphere  adequately  represented  the  mean  zonal  winds  and 
Newtonian  cooling  (^  21  levels  or  Ap  ^10  mb).  Insufficient  strato¬ 
spheric  resolution  resulted  in  a  poorly  defined  polar  night  jet  and 
Newtonian  cooling  profile  and  spurious  reflection  at  the  upper  boundary. 
They  also  suggested  that  misrepresentation  of  the  structure  of  stationary 
planetary  waves  can  lead  to  the  generation  of  spurious  external,  free 
waves  of  which  the  5  day  wave  (Geisler  and  Dickinson,  1976)  has  been 
found  to  be  dominant  in  some  forecast  models  (Lambert  and  Merilees,  1978). 

The  movement  of  tropospheric  weather  systems  at  middle  and 

high  latitudes  is  greatly  influenced  by  these  stationary,  ultralong 
Manuscript  submitted  February  6,  1980. 
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waves,  particularly  wave  numbers  1  and  2.  Therefore  a  good  forecast 
of  the  former  requires  a  good  forecast  of  the  latter.  These  planetary 
scale  eddies  are  formed  in  the  troposphere  through  baroclinic,  diabatic, 
and  orographic  processes.  During  winter  these  waves  propagate  ver¬ 
tically,  perturb  the  stratospheric  circulation,  transport  heat  and 
momentum  from  mid-latitudes  into  the  polar  regions,  and  generate  sudden 
stratospheric  warmings.  Thus  it  is  imperative  to  incorporate  a 
sufficiently  resolved  stratosphere  for  tropospheric  weather  prediction. 

In  this  report  we  document  an  effort  to  determine  the  influence 
of  stratospheric  conditions  which  are  representative  of  a  sudden  warm¬ 
ing  on  tropospheric  weather.  Quiroz  (1977),  McGuirk  (1978),  and  O’Neill 
and  Taylor  (1979)  showed  that  the  1976/77  sudden  warming  definitely 
affected  tropospheric  weather  patterns.  In  addition  Rosmond  (private 
communications,  1978)  has  noted  that  Navy  forecasts  have  lower  skill 
during  sudden  warmings.  The  goals  of  this  study  were  to  ascertain  the 
importance  of  stratospheric-tropospheric  dynamical  coupling  and  the  time 
scale  for  stratospheric  processes  to  affect  tropospheric  weather. 


II.  NUMERICAL  EXPERIMENTS 

The  UCLA  atmospheric  general  circulation  model  (Arakawa  and 
Lamb,  1977)  has  been  developed  as  a  forecast  model  at  the  Naval 
Environmental  Prediction  Research  Facility  (NEPRF)  (Payne,  1979)  and 
was  used  for  the  numerical  experiments  discussed  below.  Two  forecasts 
of  5-days  duration  each  were  performed  with  the  same  initialized  data  for 
January  19,  1977,  generated  from  NEPRF* s  global  model  initialization 
program.  During  that  week,  record  cold  temperatures  were  experienced 
over  most  of  the  continental  United  States  and  an  intense  sudden 
warming  had  developed  in  the  stratosphere.  The  reference  forecast 
had  the  usual  oj  *  0  upper  boundary  condition  at  50  mb.  The  other 
forecast  was  performed  with  an  imposed  vertical  velocity,  oj^,  at 
50  mb  based  on  the  observed  geopotential  heights  of  wave  number  1 
and  2  during  that  intense  sudden  stratospheric  warming. 

The  vertical  velocity  boundary  condition  has  two  wave  components, 

1  2 

u>I  3  u>i  +  Wj,  where  the  superscripts  represent  wavenumber,  and  are  de¬ 
fined  as: 


Wj  *  Ag(6)f1(<j>) 

to“  -  Ag(6)f2(<|>), 

Here  3  is  the  latitude  and  $  is  the  longitude,  g(Q)  is  defined  to  be 
zero  between  equator  and  30°N  and  to  be  a  sine  curve  between  30°N  and 
90°N  and  is  symmetric  with  respect  to  the  equator.  Also 
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■  cos(<t>  -  4>1) 


f2(4>)  *  cos  [(<(>  -  <J>2 ) / 2 3 


o  o 

where  0^  corresponds  to  120  W  and  ^  to  0  .  A  is  the  observed  ampli¬ 
tude  and  equal  to  -4.1  x  10  ^  mb/ sec,  or  -3.54  mb/day.  These  values 
are  crude  estimates  of  the  kind  of  vertical  velocities  which  may  be 
associated  with  sudden  warming  events.  The  values  were  obtained 

from  the  mechanistic  model  of  Schoeberl  and  Strobel  (1980).  and  ^ 

12  o 

are  chosen  such  that  and  ud^.  are  180  out  of  phase  with  their  respec¬ 
tive  components  of  the  observed  geopotential  height  to  insure  that  the 
vertical  propagation  of  wave  energy  is  upward.  Momentum  and  heat 
are  also  transported  through  the  model  top  by  the  imposed  vertical 
velocity. 


III.  NUMERICAL  RESULTS 


From  the  detailed  results  of  these  two  integrations,  we  have 
selected  three  grid  points  to  illustrate  the  salient  features.  The 
general  weather  pattern  is  shown  in  Fig.  1  which  is  day  5  of  the  inte¬ 
gration  with  imposed  upper  boundary  vertical  velocity.  Point  A  is  the 
location  of  the  maximum  negative  value  of  (=6  mb  day  at  the 
upper  boundary,  whereas  point  C  is  the  maximum  positive  value  of 
ujj  (*7.7  mb  day  ^).  Point  B,  which  coincides  with  a  deep  low  in  the 
Gulf  of  Alaska,  has  ~  0  (actually  -0.5  mb  day  ^).  To  emphasize 
stratospheric  influence  deep  in  the  troposphere,  sea  level  pressure 
differences  from  1000  mb  for  the  two  experiments  are  shown  in  Fig.  2  to 
4  for  a  5  day  period.  Perceptible  differences  are  evident  after  only 
1  day  at  points  A  and  C.  At  point  A,  after  4  days,  the  pressure 
difference  between  the  two  runs  reaches  3  mb.  At  point  B,  a  large 
change  between  day  3  and  4  reflects  the  movements  and  deepening  of 
the  low  the  pressure  difference.  A  maximum  difference  of  ~  6  mb  is 
reached  after  3.5  days  at  point  C. 
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IV.  ANALYTIC  STUDY  OF  THE  ATMOSPHERIC  RESPONSE  TO  THE 

IMPOSED  VERTICAL  VELOCITY. 

Before  interpreting  the  meaning  of  these  results,  it  is  instruc¬ 
tive  to  investigate  with  an  analytic  model  the  expected  planetary  scale 
response  of  the  atmosphere  to  the  imposed  vertical  velocity.  We  hypoth¬ 
esize  that  the  essential  features  of  the  response  can  be  model¬ 
ed  with  the  linearized  shallow  water  equations  with  a  mean  depth,  H, 
equal  to  the  equivalent  depth  of  an  isothermal  atmosphere  (~10  km): 

fkxV  +  V<j>'  *  0,  (1) 

<■*» 

gHV.V  =  Q,  (2) 

where  “  gW  =  g(h-H),  h  is  the  height  of  the  water  surface  and  is 
analogous  to  the  sea  level  pressure,  V  is  the  horizontal  velocity 
and  is  constant  with  height,  and  CQ/g)  is  the  mass  source  and  represents 
the  effect  of  imposed  vertical  velocity,  ou^. 

From  (1)  we  obtain  the  divergence  and  vorticity  equations: 


at 


at 


7.  V  -  fk  .  VxV  +  V V  =  0, 


-5-  k  .  VxV  +  fv  .  V  =  0. 
ot 


(3) 

(4) 
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Substitution  of  (4)  and  (2)  into  (3)  gives 

ot 


■^-y  V.V  +  f2V  .V  -gHV2(V  .  V)  +  V2Q  =  0. 

at 


If  0  has  wave  numbers  k  and  l  in  x  and  y  directions,  respectively,  i.e., 
Q  *  Q  exp  [i(kx  +  fcy)], 

the  Laplacian  can  be  written  as 

2  2  2  2 
V  -  -(kZ  +  JT)  h  -K  . 

If  we  use  the  notation 

D  s  V  •  V  -  De^^x+^^,  Eq.  (5)  becomes  (6) 


•^-y  D  +  (gHK2  +  f 2)D  =  K20. 

at 

For  the  initial  conditions  we  choose  V  =  <j/  =  0  at  t  =  0,  which  are 

consistent  with  the  data  initialization  for  the  numerical  experiments; 
thus. 


D  ®  0  (since  V  =  0  at  t  =  0), 


-  0  (from  (3)). 


The  solution  is 


where 


D  *  Acos(vt)  -  A, 


n  v2 

2  2  1/2 
v  -  (gHiT  +  C)LJ  . 
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Eq.  (2)  can  now  be  written 


+  gH[Acos(vt)  -  A]  =  0  , 

=  0  at  t  =  0. 

The  solution  for  <(/  is 

sln(vt)  ♦  (8HA^)tlai(k,<+tS’) 

.  fsizS  sln(vt) + 4  ^ 

L  V  V 

Since  the  surface  height  variation  is  analogous  to  the  sea  level 
pressure  variation  in  the  atmospheric  model  and  0  represents  the  im¬ 
posed  vertical  velocity  at  model  top,  Eq,  (9)  suggests  the  sea  level 
pressure  has  a  standing  oscillation  with  frequency,  v,  and  a 
linear  variation  with  time,  which  vanishes  as  f  =  0.  The  standing 
oscillation  is  a  direct  consequence  of  the  D  =  0  initial  condition 
and  the  imposed  vertical  velocity  (Q  ^  0) .  If  we  choose  D  =  -A  in¬ 
stead  for  the  initial  condition  along  with  *  0,  then  Eq.  (8)  is 
just  D  *  -A  and  the  standing  oscillation  term  in  Eq.  (9)  vanishes. 
Thus  the  mismatch  of  initial  and  upper  boundary  conditions  generates 
the  oscillatory  term  which  corresponds  to  the  creation  of  two  free 
modes  traveling  in  opposite  directions. 

A  numerical  estimate  for  is  given  below: 


Jgi(kx+2.y)  (9) 
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g  =  9.8  m/s2,  H  =  104m,  f  =  2CSin  9 =  1.3  *10~4s_1, 

P  *  7.3*10  5s  1,  0  =  a=6.4*10^m,  earth's  radius, 

.  2rr  n  2n  t  -  n  T  2rra 

k  =  —  ,  £  =  —  ,  L  =  2rracos0  =  rra,  L  =  — —  , 

L  LX  V  3 

X  V 


Thus  , 


2  3 

k  s  —  <?  =  — 

a  ’  a’ 


K2  =  k2  +lZ  *  =  3.2  x  10"13m  2, 

v  =  (gHK2+f2)1/2  =  2. 2xl0-4s”\ 


2  2 

=  3.1  x  103s  ,  Ar  =  0.32,  T  =  —  =  2.9  x  104  s  , 
i  2  v 

v  v 


Q  *  maxjuu^l  =  4.  lx  10  ^  mb  s  ^  for  each  wave  number. 


Thus  the  amplitude  of  the  standing  oscillation  is  0.17  mb  and 
its  period  is  8  hr.  The  growth  rate  of  the  second  term  in  Eq .  (9) 
is  1.2  mb  day  ^ . 

Fig.  5  shows  the  sea  level  pressure  difference  at  three  locations 
in  the  first  one  and  half  days  between  the  two  runs,  Pu  -  Pr ,  where 
Pu  is  the  sea  level  pressure  of  the  run  with  imposed  upper  boundary 
vertical  velocity,  and  Pr ,  that  of  the  reference  run.  Since  Pu  -  Pr 
is  analogous  to  the  <(/  in  the  shallow  water  equation,  the  time 
variation  of  Pu  -  Pr,  shown  in  the  figure  as  an  oscillation  super¬ 
imposed  on  a  linear  change,  corresponds  well  with  that  predicted  by 
the  shallow  water  model. 
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V. 


DISCUSSION 


The  study  shows  that  the  difference  between  the  two  simulations 
can  be  qualitatively  explained  by  Eq.  (9).  In  the  shallow  water 
equation  case,  the  oscillatory  part  of  the  solution  in  Eq.  (9)  can 
be  eliminated  by  using  a  different  initial  condition,  but  not  the 
linear  term. 

However  in  a  model  with  topographic  forcing  and  differential 
heating  due  to  sea-land  contrast  the  imposed  vertical  velocity  can 
be  initialized  as  part  of  the  stationary  waves  and  then  both  terms 
in  Eq.  (9)  are  eliminated.  In  future  proposed  experiments  the  imposed 
vertical  velocity  upper  boundary  condition  will  be  initialized  as 
part  of  the  stationary  waves  because  it  is  kept  constant  in  time.  It 
is  anticipated  that  the  difference  in  sea  level  pressure  between  such 
a  proposed  experiment  and  the  reference  experiment  will  not  be  domi¬ 
nated  by  the  direct  effect  of  the  imposed  vertical  velocity  and  will 
probably  be  smaller  in  amplitude  than  the  results  presented  in  Fig. 

2-4.  However  the  present  experiment  does  demonstrate  that  an  externally 
generated  stratospheric  disturbance  as  large  as  a  sudden  warming  can 
generate  an  observable  response  in  sea  level  pressure  with  a  time 
constant  of  the  order  of  a  day. 


I 


Acknowl ed  gmen t  s 

This  work  was  supported  by  the  NAVAIR  block  fund  under 
subtask  number  9F52-551-792 . 


REFERENCES 


Arakawa,  A*  and  V.  R.  Lamb,  1977:  Computational  design  of  the  basic 
dynamical  processes  of  the  UCLA  general  circulation  model. 

Methods  in  Computational  Physics,  Vol.  17,  Academic  Press, 

New  York,  pp.  174-207. 

Chamey,  J.  G.  and  P.  G.  Drazin,  1961:  Propagation  of  planetary 
scale  disturbances  from  the  lower  into  the  upper  atmosphere. 

^J.  Geo £hvs.  Res. ,  66 ,  83-109. 

Geisler,  J.  E.  and  R.  E.  Dickinson,  1976:  The  five-day  wave  on  a 
sphere  with  realistic  zonal  winds.  *T.  Atmos.  Sci. ,  33,  63-641. 

Kirkwood,  E.  and  J.  Derome,  1977:  Some  effects  of  the  upper  boundary 
condition  and  vertical  resolution  on  modeling  forced  stationary 
planetary  waves.  Mon.  Wea.  Rev.  105 ,  1239-1251. 

Lambert,  S.  J.  and  P.  E.  Merilees,  1978:  A  study  of  planetary  wave 
errors  in  a  spectral  numerical  weather  prediction  model* 
Atmosphere-Ocean ,  16  t  197-211. 

Lindzen,  R.  S.,  E.  S.  Batten  and  J.  W.  Kim,  1968:  Oscillations  in 
atmospheres  with  tops.  Mon.  Wea.  Rev. ,  96,  133-140. 

McGuirk,  J.  P.,  1978:  Planetary  scale  forcing  of  January  1977  weather. 
Science,  199,  293-295. 

O'Neill,  A.  and  B.  F.  Taylor,  1979:  A  study  of  the  major  stratospheric 
warming  of  1976/77.  Quart .  J.  R.  Met.  Soc . ,  105,  pp.  71-92. 


12 


Payne,  S.  W. ,  1979:  A  verification  case  study  of  the  NEPRF/UCLA  global 
forecast  system.  Fourth  Conf.  on  Numerical  Weather  Prediction, 

AMS,  Silver  Spring,  Md.  Oct.  29-Nov.  1,  1979,  155-160. 

Quiroz,  R.  S. :  1977,  The  tropospheric-stratospheric  polar  vortex 
breakdown  of  January  1977.  Geophvs .  Res.  Lett . ,  4_,  151-154. 

Schoeberl,  M.  R.  and  M.  A.  Geller,  1977:  A  calculation  of  the  structure 
of  stationary  planetary  waves  in  winter.  J_.  Atmos.  Sci. ,  34 , 
1235-1255. 

Schoeberl,  M.  R.  and  D.  F.  Strobel:  1978a,  The  zonally  averaged  cir¬ 
culation  of  the  middle  atmosphere.  J.  Atmos .  Sci. .  35 .  577-591. 


Sea  level  pressure  -  1000  mb  at  point  A 


-0.2  mb/hr 


►  WITH  IMPOSED  o> 
•CONTROL 


-30.9 


A  - 


32.2 


Fig.  3  -  Sea  level  pressure  -  1000  mb  at  point  B 
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Fig.  5  -  Detailed  sea  level  pressure  difference  between  the 
two  runs  at  point  A,  3,  and  C  for  the  first  1.3  days 
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